Coherent exciton transport and trapping on long-range interacting cycles 
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We consider coherent exciton transport modeled by continuous-time quantum walks (CTQWs) 
on long-range interacting cycles (LRICs) , which are constructed by connecting all the two nodes ol 
distance m in the cycle graph. LRIC has a symmetric structure and can be regarded as the extensions 
of the cycle graph (nearest-neighboring lattice). For small values of m, the classical and quantum 



return probabilities show power law behavior p(t) 



t~ - 5 and ir(t) 



respectively. However, 



for large values of m, the classical and quantum efficiency scales as p(t) ~ t~ and n(t) ~ t~ . We 
give a theoretical explanation of this transition using the method of stationary phase approximation 
(SPA). In the long time limit, depending on the network size N and parameter m, the limiting 
probability distributions of quantum transport show various patterns. When the network size N 
is an even number, we find an asymmetric transition probability of quantum transport between 
the initial node and its opposite node. This asymmetry depends on the precise values of N and 
m. Finally, we study the transport processes in the presence of traps and find that the survival 
probability decays faster on networks of large m. 

PACS numbers: 05.60.Gg, 03.67.-a, 05.40.-a 
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I. INTRODUCTION 

In the past few years there has been a growing inter- 
est in continuous-time random walks (CTRWs) [3, H, E| . 
The particular surge of increasing interest can be partly 
attributed to its close connection with the classical dif- 
fusion modeled by the tight-binding model in condensed 
matter [J] . The quantum mechanical analog of the clas- 
sical diffusion process defined on complex networks has 
also been studied with respect to the localization der- 
ealization transition in the presence of site disorder ||. 
In the literature, there are two main types of quan- 
tum walks: continuous-time and discrete-time quantum 
walks [6] . Discrete-time quantum walks evolve by the ap- 
plication of a unitary evolution operator at discrete time 
intervals, and continuous-time walks evolve under a time- 
independent Hamiltonian Q ■ It has been shown that on 
some graphs, propagation between two properly chosen 
nodes is exponentially faster in the quantum case Q. In 
this respect, quantum walks provide a good framework 
for the design of quantum algorithms in the application 
of quantum computation . 

Here, we focus on continuous-time quantum walks 
(CTQWs) . Previous work have studied CTQWs on some 
particular graphs, such as, the line cycle hy- 

percube [13], Cayley tree [H, [3], dendrime rs fl5j and 
other regular networks with simple topology [lg, [lZ| • In 
Ref. [l8j . the authors studied the coherent exciton dy- 
namics on discrete rings under long-range step lengths 
distributed according to i? -7 (7 ^ 2). The strength of 
the long-range interaction is a power law decay of the 
distance of the nearest-neighboring lattice. They find 
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that the long-range interactions give no influence to the 
efficiency of the coherent exciton transport [l8[ . 

In this paper, we study the effect of long-range in- 
teractions on a new network model, namely long-range 
interacting cycles (LRICs). LRICs are constructed by 
connecting all the two nodes of distance m in the cy- 
cle graph (nearest-neighboring lattice). Therefore, the 
network model has a symmetric structure and can be 
regarded as the extensions of the cycle graph (nearest- 
neighboring lattice). The newly added edges with large 
m are long-range interactions and serves as shortcuts in 
the nearest-neighboring cycle graph. A detailed descrip- 
tion of the network structure will be given in the next 
section. 

Since the structure of LRICs is completely symmetri- 
cal as the nearest-neighboring cycle graph, we are able to 
analytically predict the dynamical behavior of the coher- 
ent and incoherent transport. The paper is organized as 
follows: In Sec. II we give a description to the structure 
of LRICs. In Sec. Ill, we briefly review the properties of 
CTQWs on general graphs. In Sec. IV, we derive ana- 
lytical results for LRICs and study the efficiency of the 
classical and quantum transport by considering the scal- 
ing of the return probability. Long time averages of the 
transition probabilities are also studied in this section. 
In Sec. V, we study trapping process on LRICs. Con- 
clusions and discussions are given in the last part, Sec. 
VI. 



II. TOPOLOGY AND STRUCTURE OF LRICS 

Long-range interacting cycles (LRICs) can be con- 
structed as follows: First, we construct a cycle graph 
of N nodes where each node connected to its two near- 
est neighbor nodes. Second, two nodes of distance m 
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(a) N=10 (b) N=10 

m=2 m=3 

FIG. 1: (Color online) Long-range interacting cycles G(10, 2) 
(a) and G(10, 3) (b). 

in the cycle graph are connected by additional bonds. 
We continue the second step until all the two nodes of 
distance to have been connected. Hence, the LRICs, 
denoted by G(N,m), are characterized by the network 
size N and long-range interaction parameter m. LRIC is 
a one-dimensional lattice with periodic boundary condi- 
tions and all nodes of the networks have four bonds. The 
structure of G(10, 2) and G(10, 3) is illustrated in Fig. □ 
It is interesting to note that all the LRICs have the 
same value of connectivity k = 4 and the parameter to 
adjusts the interaction range of the cycles, thus LRICs 
provide a good facility to study the effects of long-range 
interaction on the transport dynamics. 



equation [lj| [lj| [2(|. Then the classical and quantum 
transition probabilities to go from the state \j) at time 
to the state \k) at time t are given by pt i(t) = (fc|e~ tA |.7') 
and w hj (t) = \a kj {t)\ 2 = |(fc| e -^|j)|2 [HG3, res P ec - 
tively. Using E n and \q n ) to represent the nth eigenvalue 
and orthonormalized eigenvector of H , the classical and 
quantum transition probabilities between two nodes can 
be written as [H Qj, [H, Hj 

PkAt) = T / e ~ tEn ( k M(ln\j), (1) 



7r*j(t) = \a k At)\ 2 = \E n e- UE -(k\q n )(q n \j)\ 2 f 
= ^ nl e-^- E ')(k\ qn )( qn \j)(j\ qi )( qi \k). {Z > 

For finite networks, 7Tfcj(i) do not decay ad infini- 
tum but at some time fluctuates about a constant value. 
This value is determined by the long time average of 

Xk,j = limr^oo y J T TTfcj (t)dt 

= T, n ,i( k \<ln)(qn\j)(j\qi)(qi\k) 

= J2nA S E n ,E l (k\q n ){qn\j){j\qi)(qi\k). 

where &E n ,E l takes value 1 if E n equals to Ei and oth- 
erwise. Generally, to calculate pk,j(t), nk,j(t) and Xk,j 
all the eigenvalues E n and eigenvectors \q n ) are required. 
For some regular graphs, the eigenvalues and eigenvectors 
can be analytically obtained. In the following section, we 
find analytical results of the eigenvalues and eigenstates 
for LRICs, and calculate these quantities according to 
the above Equations. 



III. COHERENT EXCITON TRANSPORT ON 
GENERAL GRAPHS 

The coherent exciton transport on a connected net- 
work is modeled by the continuous-time quantum walks 
(CTQWs), which is obtained by replacing the Hamilto- 
nian of th e sy stem by the classical transfer matrix, i.e., 
H = -T [Il[l9||. The transfer matrix T relates to the 
Laplace matrix by T — ~"fA, where for simplicity we as- 
sume the transmission rates 7 of all bonds to be equal and 
set 7 = 1 in the following [l^, [lU H3| • The Laplace ma- 
trix A has nondiagonal elements Aij equal to — 1 if nodes 
i and j are connected and otherwise. The diagonal el- 
ements An equal to degree of node i, i.e., An — ki. The 
states \j) endowed with the node j of the network form 
a complete, ortho-normalised basis set, which span the 
whole accessible Hilbert space. The time evolution of a 
state \ j) starting at time to is given by \ j, t) = U(t, 
where U(t,to) — exp[~iH{t — t^)] is the quantum me- 
chanical time evolution operator. The transition ampli- 
tude akj(i) from state \j) at time to state \k) at time 
t reads a>k,j(t) — (k\U(t,0)\j) and obeys Schrodingers 



IV. COHERENT TRANSPORT ON LRICS 

A. Analytical results 

In the subsequent calculation, we restrict our attention 
on the graph of long-range interacting cycles (LRICs). 
The network organizes in a very regular manner and has 
a periodic boundary condition. The Hamiltonian matrix 
H of G(N, m) (to 2) takes the following form, 

r 4, if i = j, 

l. 0, Otherwise. 

And the Hamiltonian acting on the state \ j) can be writ- 
ten as 

H\j) = 4|j) - \j - 1) - |j + 1) - |j - to) - \j + to). (5) 

The above Equation is the discrete version of the Hamil- 
tonian for a free particle moving on the cycles. Using the 
Bloch function approach for the periodic system in solid 
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state physics (2l|, the time independent Schrodingers 
Equation reads 



H\^ n )=E n \^ n 



(6) 



The Bloch states \ipn) can be expanded as a linear com- 
bination of the states \j) localized at node j, 



\1>n) 



1 N 

3 = 1 



ii)- 



(7) 



Substituting Eqs. (O and (0 into Eq. ([6]), we obtain the 
eigenvalues (or energy) of the system, 



E n = 4 — 2 cos 9 n — 2 cos(m6 n ) 



(8) 



The periodic boundary condition for the network requires 
that the projection of the Bloch state on the state |iV+l) 
equals to that on the state |1), thus 8 n = 2nir/N with n 
integer and n € [1, N}. Replacing \q n ) by the Bloch states 
\xj) n ) in Eqs. (JTJ) , J2J and (|3|), we can get the classical and 
quantum transition probability 



e -tE ne -i(k-j)2nw/N 



ilE„ 



-i(k-j)2mr/N\2 



~~ I AT En 1 

and the long time averages of ^kj(t) is given by 

-i(fc-j)(n-i) 27I 7 iv 



Xfc, 



E n ,E,e 



(9) 



(10) 



(11) 



>i. i 



Interestingly, when k = j, the transition probability 
is reduced to the return probability, which means the 
probability of finding the exciton at the initial node. In 
Ref. [22| , the authors use the return probability to quan- 
tify the efficiency of the transport. In the next subsec- 
tion, we will analyze return probability and try to com- 
pare the efficiency between the classical and quantum 
transport. For our regular cycles, the return probability 
is independent on the initial node. The average return 
probability can be written as, 



Pit) 



-tEn 



and 



<t) =^E^-(*) = Ki(*)l a = |a(*)l s 

= l££„e-^| 2 . 



(12) 



(13) 



Eqs. (1121) and (Ti~3|) hold for finite networks. For infinite 
networks, i.e., N — ► oo, the 9 values are quasi-continuous 
in Eq. ([5]). In the continuum limit, on one hand, the 
eigenvalues of Eq. ([8]) can be rewritten as, 



E m {9) = 4- 2 cos 9 - 2 cos m9. 



(14) 




FIG. 2: (Color online) Eigenvalues E m (8) versus 8 for m — 2 
(a) and m = 10 (b). For m = 2, there are two maxima 
and minima in [0, 2ir), which are indicated by the arrows in 
the plot (Fig. 2 (a)). For m — 10, there are ten maxima 
and minima in [0, 27r), and the maximal and minimal points 
are indicated as E max (8) and E min {8) (See dashed curves in 
Fig. 2 (b)). 



on the other hand, the classical and quantum return 
probabilities in Eqs. (fT2|) and (|13p can be written as the 
following integral form, 

1 f 2 " 

Pm(t) = ^- / ex P (-tE m (9))d9, (15) 
27r Jo 



and 



1 



7r w (t) = \a(t)\ 2 = I— / exp{-itE m {9))d9\ 2 . (16) 
27T Jo 

Fig. [2] shows E m (9) versus 9 for m = 2 (a) and m — 
10 (b). We note that E m (9) is an oscillatory function, 
and there are more regular oscillations for large values of 
to. The number of maxima (or minima) of E m {9) in the 
range [0, 27r) is m. As we will show, these extreme points 
give contributions to the integrals when we calculate the 
classical and quantum efficiency in Eqs. (fT5|) and (fTB) . 



B. 



Efficiency and scaling of the classical and 
quantum transport 



In this subsection, we consider the efficiency of the 
classical and quantum transport. We calculate the in- 
tegrals of Eqs. (jTSJ) and (flB)) using the stationary phase 
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FIG. 3: (Color online) Return probability p m (t) and n m (t) for 
different values of m. (a) Classical return probabilities p m (t) 
for m — 2, m = 10 and m = 100. The black curves are ex- 
act results obtained according to Eq. (|12[) . green curves are 
the analytical results obtained using stationary phase approx- 
imation (SPA) (See Eqs. lfB4)l and (fB6)l in the Appendix), 
(b) Quantum return probability ^(t). The solid curve is the 
exact result and the dashed curve is the analytical result ob- 
tained using stationary phase approximation (See Eq. (|C1|I in 
the Appendix). (c)7rio(t) versus t. The solid curve is the ex- 
act result and the dotted line is the power law t~ . (d)7Tioo(t) 
versus t. The solid curve is the exact result and the dashed 
curve is analytical prediction according to Eq. (|19|l . All the 
exact results are obtained from LRICs of size N = 10000. 



approximation (SPA) (See Appendix A). We find that 
the classical and quantum return probabilities show dif- 
ferent scaling behavior for small values and large values 
of m. 

For small values of to, we get an asymptotical expres- 
sion for the classical p m (t), 

p m (t)n—}-=~t-™. (17) 
Zmyirt 

(See derivation in Appendix B). For large values of m, 



we also get an approximate result (See Appendix B), 

Mt)*^-* -1 . (18) 

We note that p m (t) scales as t~° for small values of to, 
however, for large to, the scaling becomes as p m (t) ~ t . 

Quantum mechanically, we also find different scaling 
behavior of ir m {t) for small value and large value of to. 
For small values of m, 7r m (t) is an oscillatory function 
multiplied by X/t (See Eq. (|C1[) in Appendix C for to = 
2). For large values of to, there is an approximate result 
given by Eq. (|C8[) in Appendix C, 

. , sin 2 At 9 . „. 

4^5"-* • ( 19 ) 

Therefore, quantum transport of small m displays the 
same scaling behavior n(t) ~ t^ 1 while transport of large 
to shows scaling w(t) ~ t~ 2 . It is interesting to note 
that, both for the classical and quantum transport, the 
scaling exponent for large to is twice the exponent for 
small values of to. This is one of the main conclusions in 
this paper. 

In order to test the theoretical predictions, Fig.[3]shows 
the classical and quantum return probabilities for LRICs 
of N = 10000 with to = 2, to == 10 and 100. Fig. [3] 
(a) shows the classical return probability. We note that 
P2{t) and Pxo(t) displays the same scaling t -0 5 , but for 
to = 100 the scaling becomes pioo(t) ~ t~ X - The results 
are in good agreement with the analytical predictions of 
Eqs. (fTT)) and (fP8"|) . Fig. [3] (b) shows the quantum n m (t) 
for to = 2 and the analytical prediction of Eq. (|C1[) in 
Appendix C. Both the results exhibit power law ^(i) ~ 
i -1 . The same scaling behavior (< _1 ) is also observed for 
to = 10 (See Fig. (c)). In Fig. (d), we show ir 100 (t) 
and the analytical result predicted by Eq. (JTUJ) . Both the 
results display the same scaling t~ 2 . 



C. Long time averages on finite networks 

In this section, we consider the long time averaged 
transition probabilities on finite networks. Classically, 
the long time liming probabilities equal to the equip- 
partitioned probability 1/N [23]. Quantum mechanically, 
the long time averages of the transition probabilities does 
not lead to equip-partition. For LRICs, the long-time 
averaged probability is determined by Eq. (jlip but the 
distribution patterns are complex for different network 
parameters N and to. For the cycle graph (nearest neigh- 
boring lattice), the limiting probability distribution de- 
pends on the parity of the network size N. Fig. 2] shows 
the distribution patterns of the limiting transition prob- 
ability on networks of N = 100 with various values of 
in. The initial excitation is located at node 1. As we 
can see, there are high probabilities to find the exciton 
at the initial node 1 and the opposite node 51, this fea- 
ture is a natural consequence of the periodic boundary 
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FIG. 4: Long-time averaged probability distribution \k, 1 for 
CTQWs on networks of size N = 100 with different values of 
m. 



FIG. 5: Quantum mechanical limiting probabilities \k, 1 on 
networks of size N = 75 with different values of m. 



condition of the graphs [20| . For odd- numbered networks 
TV e Odds, there is a higher probability to find the initial 
node than that at other nodes [2(J. Fig. [5] shows the dis- 
tribution patterns for networks of N = 75 with various 
to. The patterns depend on the specific network param- 
eters and there are high probabilities to find the exciton 
at some particular nodes. We also note that the patterns 
of Xk,i are the same for some different values of to, this 
feature can be explained by the identical degeneracy dis- 
tribution of the eigenvalues for different values of to [2(| • 

It is worth mentioning that for even-numbered net- 
works, there are high probabilities to find the exciton 
at the initial node and opposite node. For networks of 
N = 100, we find the two probabilities are exactly equal 
to each other for all the values of to, i.e., Xi,i — X5i,i- 
However, for some other even-numbered network size N, 
this is not true (20| . For some particular values of N 
and m, for instance N — 108 and to = 2, the probabil- 
ity of finding the exciton at the initial node differs from 
the probability of the finding the exciton at the opposite 
node. Such asymmetry is small and not easy to be ob- 




FIG. 6: (Color online) A(l,JV/2) = ( X i,i - Xiv/2+i,i)/(xi,i + 
Xiv/2+1.1) as a function of the network size iV for different 
values of m. The solid line indicates the power law decay 
A(l, N/2) ~ TV -1 . 
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served from the limiting probability distributions [2CJ • To 
detect such asymmetry of the probabilities, we define the 
quantity A (l,N/ 2) = (xi,i - Xw/2+i,i)/(xx,i+Xjv/2+i,i) 
as a function of the network size N for different values 
of to |20| | . The asymmetry is indicated by the nonzero of 
this quantity while A(l,iV/2) = corresponds to iden- 
tical values of xi,i an d Xn/2+i,i)- A plot of A(l,N/2) 
versus N for to = 2, to = 3 and to = 4 are shown in 
Fig. [6] We find that the points break into several clus- 
ters, whereas some clusters A(l, N/2) decreases with the 
network size N as a power law: A(l, N/2) ~ iV _1 [20| . 

As we have shown in Fig. (6[ the asymmetry appears 
at some particular values of N and to. However, we are 
unable to predict which particular parameters of N and 
m are related to such asymmetry. This is an interesting 
issue and similar phenomena is also found in Ref. [20j. 



V. TRAPPING ON LRICS 

An im por tant process related to random walk is trap- 
ping [24I |25| . Trapping problems have been widely stud- 
ied in the frame of physical chemistry, as part of the gen- 
eral reaction-diffusion scheme (26[. Previous work has 
been devoted to the trapping problem on discrete-time 
random walks [U H^|- However, even in its simplest 
form, trapping was shown to yield a rich diversity of re- 
sults, with varying behavior over different geometries, di- 
mension, and time regimes [28j . The main physical quan- 
tity related to trapping process is the survival probabil- 
ity, which denotes the probability that a particle survives 
during the walk in a space with traps. 

In this paper, we consider trapping using the approach 
based on time dependent perturbation theory and adopt 
the methodology proposed in Ref. [29| . In Ref. [2^|, the 
authors consider a system of N nodes and among them 
M are traps (M < N). The trapped nodes are denoted 
them by m, so that to e M. The new Hamiltonian 
of the system is H = Ho + iT, where Ho is the orig- 
inal Hamiltonian without traps and iT is the trapping 
operator. T has to purely imaginary diagonal elements 
T mm at the trap nodes and assumed to be equal for all 
to (r mm = r > 0). See Ref. d| for details. The new 
Hamiltonian is non-hermitian and has N complex eigen- 
values and eigenstates {Ei, \^i}} (I — 1,2, ...,N). Then 
the quantum transition probability is 

1 

where (^;| (I — l,2,...,N) is the conjugate eigenstates 
of the new Hamiltonian. In order to calculate iTkj(t), 
all the complex eigenvalues and eigenstates {Ei, \SSfi)} 
(I = 1,2,...,N) are required. Here, we numerically cal- 
culate TTk,j (t) by diagonalizing the Hamiltonian H using 
the standard software package Mathematica 5.0. 

Equation (|2"D|) depends on the initially excited node j. 
The average survival probability over all initial nodes j 




FIG. 7: (Color online) Survival probabilities IIm(£) and Pm(£) 
for LRICs of JV = 100 and different values of m. In the 
calculation, five trap nodes are randomly selected from the 
cycles and we set F = 1. The curves are averaged over distinct 
trapping realizations. 

and all final nodes k, neither of them being a trap node, 
is given by, 

HmW = iv^i? £ £ (21) 

For continuous-time random walks (CTRWs), we in- 
duce trapping analogously as the CTQWs, where the 
new transfer matrix is modified by the trapping matrix 
as T = To — r. The mean survival probability analogous 

to Eq. dsn is p M (t) = j±m E^m YtkiMPki® m- 

Fig. [7] shows the quantum and classical survival prob- 
abilities on LRICs of N = 100 with m = 2, to = 5 
and to = 10. Five trapped nodes are randomly selected 
from all the nodes and T = 1 is fixed in the numerical 
calculation. For each specific trapping configuration, we 
calculate the survival probability and average it over dif- 
ferent configurations. As we can see from Fig. [71 both 
the quantum and classical survival probabilities decays 
fast on LRICs with large values of m. This is opposite to 
the case in Ref. [30( where long-range interaction leads 
to a slower trapping of the excitation. 

VI. CONCLUSIONS AND DISCUSSIONS 

We have studied coherent exciton dynamics modeled 
by continuous-time quantum walks (CTQWs) on long- 
range interacting cycles (LRICs). We have shown that 
both the efficiency of the classical and quantum trans- 
port display power laws, and the exponents for LRICs 
with large values of to are twice the exponents of LRICs 
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with small values of m. Theoretical calculation of the 
return probability using stationary phase approximation 
supports this finding. In the long time limit, the limiting 
probability distributions of quantum transport show var- 
ious patterns on finite size networks. When the network 
size N is an even number, we find an asymmetric transi- 
tion probability between the initial node and its opposite 
node. This asymmetry depends on the precise values of 
N and to. Finally, we study trapping process on LRICs 
and find that long-range interaction (large m) leads to a 
fast decay of the survival probability. 

It is worth mentioning that the return probability dis- 
plays different scaling behavior for small values and large 
values of m. However, we did not give a quantitative re- 
lation between the scaling exponent and the parameter 
to. We only know the scaling behavior for the thresh- 
olds of small and large to, the scaling behavior in the 
medial region of m is still unknown. In addition, the 
limiting probability distributions show various patterns 
on finite size networks. These patterns is a natural re- 
sult of the interference phenomena in coherent transport 
on finite systems. The asymmetry of the limiting prob- 
ability on even-numbered networks is also an interesting 
and strange feature of quantum walks, which deserves 
our further investigation [23|, [3l| . The long-range inter- 
action in LRICs leads to a fast exciton trapping and is 
opposite to the conclusions in Ref. [30| . We also note 
that the quantum return probability in Ref. [3| scales 
the same behavior w(t) ~ t^ 1 for various long-range in- 
teractions (7 2). The different behavior of trapping 
and transport efficiency may be caused by the distinct 
type of long-range interactions. 
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The SPA asserts that the main contribution to this in- 
tegral comes from those points where f{x) is stationary 
[df(x)/dx = f'(x) = 0]. If there is only one point Xq for 
which f(x ) = and d 2 f(x) /dx 2 \x = f"(x ) ^ 0, the 
integral is approximated asymptotically by, 



7(A) 



1 



-A/(so) 



(A2) 



If there are more than one stationary points satisfy 
[df(x)/dx = f'(x) = 0], then the integral /(A) is ap- 
proximately given by the sum of the contributions [each 
being of the form given in Eq. (|A2j) ] of all the stationary 
points [HI. 



APPENDIX B: CALCULATION OF THE 
CLASSICAL p m (t) USING SPA 

We apply SPA to calculate the classical p m (t) in 
Eq. (fT5|) . The stationary points of this integral sat- 
isfy E' m (9) — 2sin# + 2m sin m8 — 0. The number of 
stationary points equals to 2m (m maxima and m min- 
ima, see Fig. [5]) in the range 9 £ [0, 27r). We denote m 
maxima stationary points as 9™ ax (i = 1, 2, m) and m 
minima stationary points as 9™" m (% = 1, 2, ...,m). Then 
the integral of Eq. (jTSj) yields, 



n (t) ps V m 1 -tE^iev 1 " 1 ) 

1 „-t£ m (er ax ) ^ ' 



^/27rt£^(e™' lx ) 



which is mainly determined by the small values of E m (9i). 
Considering E m (0 min ) < 4 < E m (9 max ), contributions 
from the maximal stationary points in the above equation 
is negligible. Therefore, p m (t) can be simplified as, 
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APPENDIX A: THE STATIONARY PHASE 
APPROXIMATION (SPA) 

Stationary phase approximation (SPA) is an approach 
for solving integrals analytically by evaluating the inte- 
grands in regions where they contribute the most [ToL [l8l . 
|32| . This method is specifically directed to evaluating 
oscillatory integrands, where the phase function of the 
integrand is multiplied by a relatively high value. Sup- 
pose we want to evaluate the behavior of function /(A) 
for large A, 

/(A) = i- J e- Xf{x) dx. (Al) 



>(*) - E 



1 



(B2) 



For small values of m, the global minimum E m {&) at 
9 = is sufficiently separated from (smaller than) other 
local minima. The sum in Eq. (|B2[) is mainly from the 
contribution at the global minimum 9 = 0, thus, 



Pm(t) 



1 



-tE m (Q) 



(B3) 



Substituting the relation E" n (0)|e=o = 2 + 2m 2 and 
E m (0) = into Eq. (|B3]). we get, 



Pm(t) 



\J\-Kt{\ + m 2 ) 2mV7rt 



(B4) 



For large values of m, the global minimum E m (9) at 
9 = is not sufficiently separated from (smaller than) 
other local minima. The sum in Eq. (|B2jl contains contri- 
butions from all the minimal stationary points. Noting 
that for large values of to, the stationary points 9™" ax 



8 



and QV n - m are approximately equidistant, i.e., #™ m w 
2(i - l)n/m, 9™ ax w (2i - 1)tt/to (i = l,2,...,m). 
Therefore, we get the approximations: E' 1 ! n (6 mm ) ss 2m 2 , 
E'^(9 max ) w -2m 2 . Thus Eq. jB2| can be written as, 



Pm{t) 



1 



2aAtI 



e(t), 



(B5) 



In the contin- 



where Q(t) = ±^1 
uum limit of large m, Q(t) equals to the integral 
h Jo' T e~ tEmt " (e) dO, where E min {9) = 2 - 2cos6> (See 
the dashed curves in Pig.[5](b)). We apply the method of 
SPA again to evaluate this integral and find that the con- 
tribution is mainly from the stationary point 9 = 0, which 
lead to 6(i) w -^=. The classical p m (t) of Eq. (JBSJ 
transforms into, 

(B6) 



APPENDIX C: CALCULATION OF THE 
QUANTUM iv m (t) USING SPA 

We calculate the integral of quantum return probabil- 
ity in Eq. (fT6|) . We also find the integral displays different 
scaling behavior for small values and large values of m. 

For small values of m, we consider the case m = 
2, where there are four stationary points: 9"" n = 0, 

= 2vr - 



ymax 



ArcCos(-l/4), 9? 



ArcCos(—l/4) (See Fig. [2] (a)). The second-order deriva- 
tions at these points yield E"(9? in ) = 10, E"(9^ in ) = 6 
and E"(9f ax ) = E"(9 7 2 nax ) = -15/2. The corresponding 
spectral eigenvalues at the four points are E(9™ m ) = 0, 
E{6f in ) = 4 and E(9™ ax ) = E(9™ ax ) = 25/ '4. Using the 
method of SPA, we obtain the integral of Eq. (fTB|) for 
in = 2 as, 



7r 2 (i) =|^ F / Q ihr eipHtB TO (fl))de| s 



-it-0 



1 V2ttiM0 
1 



-it-25/4|2 



^/27rit-(-15)/2 

i 3^(12 + \/l5cos4* - 4v / 5sin9t/4 
-4V3sin25t/4) - f 1 . 



(CI) 



Therefore, the quantum mechanical efficiency scales as 
7T2(t) ~ i" 1 for m = 2. For other small values of m, 
the calculation is analogous. The result is also an oscil- 
latory function multiplied by 1/t. This suggests that the 
quantum transport of small m displays the same scaling 
behavior n(t) ~ t . 



For the case of large m, the integral of Eq. (|T6|) comes 
from 2m stationary points, 



y/2TritE' J ^(8^ in ) 



its m (e s m '") 
-itE m (0j naa ') 



^/27rit£ I '4(e™ a: 

J J_ V^™ e -itE m (e m ' n ) 

' 71 Z^i=l L 



(C2) 



-itB m (e™ ax )i2 



1 2\Arii m 

+ I J_v m e" 

2V-7rit rrt £-u=l 

where in the last approximation E'^ n (9 max ) as —2m 2 and 
E" n (9 mm ) w 2m 2 is applied. In the continuum limit of 
m — ► oo, the sum in the above equation can be written 
as the integral form, 



TT7 ^ ' 



-itE m (6 



2tt 



and 



(C3) 



-it£ m (e™ ax ) _ J_ 
2tt 



2tt 



exp(~itE max (9))d9. 



(C4) 

Noting that E mm (9) = 2 - 2cos6» and E mm {9) = 6 - 
2cos6> (See the dashed curves in Fig. H (b)), Eq. (|C2|) 
can be rewritten as, 



7T m (t) 



— i f n 2,r e:cp(-ii(2 - 2 cos 6>))d<2 

1 2v / irit 2tt JO V " (C5) 

— C exp{~it{Q^2cos9))d9\ 2 . K 1 



■nit 27T JO 




For the two integrals in the above equation, we apply 
SPA again and find that the contribution of this integral 
is mainly from two stationary points 9 = and 9 = n. 
Thus 



1 f 2 * 

— / exp(-it(2- 2 cos 9))d9 
2tt Jo 



1 



1 



2^/^^ii 2^-mt 



and 



(C6) 



1 



-Ait 



±- I esp(-tf(6-2cos0))d0i- 

^tt Jo 2\/mt 2y/—Trit 

(C7) 

Substituting these relations into the Eq. (|C5[) . we get, 



8;'/ 



TTm (t) 



sin 2 4± 



(C8) 
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